a=load('cube_MoM_HH.dat');
b=load('cube_EPA_HH.dat');
tick=[0:1:180];

figure(1);
plot(tick,a,'b-');
hold on;
plot(tick,b,'r+','MarkerSize',3);
xlabel('\theta_s');
ylabel('RCS(dBsm)');
legend('MoM','EPA');
title('HH polarization');